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A minimal, analytically manageable Galerkin type model for convection in binary mixtures subject 
to realistic boundary conditions is presented. The model elucidates and reproduces the typical 
bifurcation topology of extended stationary and oscillatory convective states seen for negative Soret 
coupling: backwards stationary and Hopf bifurcations, saddle node bifurcations to stable strongly 
nonlinear stationary and traveling wave (TW) states, and merging of the TW solution branch 
with stationary states. Also unstable standing wave solutions are obtained. A systematic analysis 
of the concentration balance for liquid mixture parameters has lead to a representation of the 
concentration field in terms of two linear and two nonlinear modes. This truncation captures the 
important large-scale effects in the laterally averaged concentration field resulting from advective 
and diffusive mixing. Also the fact that with increasing flow intensity along the TW solution branch 
the frequency decreases monotonically in the same way as the mixing increases — the variance of the 
concentration distribution decreases — is ensured and reproduced well. Universal scaling relations 
between flow intensity, frequency, and variance of the concentration distribution (degree of mixing) 
in a TW are predicted by the model and have been confirmed by numerical solutions of the full 
equations. The validity of the model is checked by comparison with numerical solutions of the full 
field equations. 

PACS number(s): 47.20.-k, 47.10. +g, 03.40. Gc 
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I. INTRODUCTION 

A great deal of effort has recently been undertaken 
to investigate convection in binary fluid mixtures as an 
example for pattern formation far from equilibrium [Q. 
This system provides an experimentally convenient de- 
vice with a well established theoretical description 



allowin, 
investigations ||15 



quantitative comparisons of theoretical 
"29! with experiments. For a review 
and additional references, see |p],p^. Compared to con- 
vection in ordinary one-component fluids the spatiotem- 
poral properties are far more complex due to the influence 
of Soret sustained concentration gradients. The evolu- 
tion of the concentration fleld is governed by the inter- 
play of typically strong nonlinear convective transport 
and mixing, weak dissipative solutal diffusion, and the 
Soret effect ^, p"3Ul4|| . The latter is a source of concen- 
tration fluctuations. It generates concentration gradients 
in response to the externally applied temperature differ- 
ence and to local temperature gradients. The strength 
of the Soret coupli ng is measured by the dimensionlcss 
separation ratio if) ||l| p3| , fl . 

The concentration field changes the convective dynam- 
ics via solutal buoyancy forces entering into the momen- 
tum balance. In this way concentration gradients directly 
influence the flow which in turn changes and mixes the 
concentration. In binary liquids, this strongly nonlin- 
ear feed back is only weakly damped by small diffusive 
homogenization so that the concentration field distribu- 
tion shows significant anharmonic and boundary layer 
structures. It is, however, the existence of the feed-back 
loop that ultimately causes convection in binary mixture 
to exhibit such a rich variety of patterns arising from 



stationary and oscillatory instabilities: Depending 
on the parameters the hydrodynamic balance equations 
show convective solutions that bifurcate out of the qui- 
escent conductive basic state in the form of (i) straight, 
stationary, parallel rolls, (ii) traveling waves (TWs) con- 
sisting of propagating rolls, (iii) standing wave (SW) os- 
cillations, and (iv) stationary squares. Besides these pri- 
mary states, there are close to onset of convection pulse 
like, spatially localized traveling wave (LTW) states con- 
sisting of only few TW rolls, oscillating square patterns 
pl| , cross-roll structures [Q, and also (spiral) defect 
chaos 

In the present paper the focus is on 2-dimensional 
(2D) spatially extended convective structures consisting 
of straight parallel rolls that occur at negative separation 
ratios —0.6 < "0 < 0. For typical fluid parameters con- 
vection arises via an oscillatory subcritical bifurcation. 
The emerging solution branch locates unstable TWs that 
are "weakly nonlinear" only near the onset. These un- 
stable waves become strongly nonlinear and anharmonic 
pO| , p^ well before the occurrence of a saddle node at 
which they are stabilized on an upper solution branch. 
Simultaneously, the TW propagation speed slows down 
from its large value at the Hopf bifurcation threshold to- 
wards zero at the final transition to steady overturning 
convection (SOC). There, the amplitude of the concen- 
tration wave vanishes since in the SOC state the fluid is 
well mixed to a mean concentration level except within 
narrow boundary layers. This SOC state is somewhat 
similar to the convective rolls in one-component fluids. 

The bifurcation topology decribed above has been ver- 
ified by several experimental groups (e. g. ||]; for addi- 
tional references see [p]JT^). A detailed insight into the 
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spatiotemporal variations of TW and SOC states along 
their upper solution branches and their parameter de- 
pendence provided numerical simulations of Barten et al. 

A quantitative description of the whole bifurca- 
tion branches ||2^,Q including the lower branches that 
were unavailable to Barten et al. was obtained recently 
with a MuLTi Mode Galerkin (MMG) expansion in- 
cluding several hundred modes. The MMG predictions 
agree very well with results from finite difference 

Marker And Cell (MAC) simulations of the full field 
equations supplemented by a control process which al- 
lows to evaluate also unstable TW and SOC states Jsst . 

Since the first observations of TWs much theoretical 
research activity has been devoted to an understand- 
ing and to develop models for these phenomena. Based 
upon an earlier few-mode Galerkin model Cross |3^] 
and also Ahlers and Liicke investigated Soret-driven 
convection with permeable boundary conditions. They 
found TW states only locally at the onset of convection 
since for permeable conditions the onset of TW convec- 
tion is tricritical. Linz et al. implemented imperme- 
able conditions yielding a subcritical primary bifurcation 
to TWs. This observation is in agreement with small 
ampHtude computations of Schopf et al. |2^. Both ap- 
proaches however do not explain the stabilization of the 
TW branch via a saddle node bifurcation. Bensimon et 
al. jl^ considered the case of weak Soret coupling by 
means of a small-V' expansion treating the concentration 
field numerically. They observed a stable TW branch 
and interpreted the TW-SOC transition as a boundary 
layer induced instability. Due to the expansion in ip and 
due to the weak diffusion limit the application range in 
fluid parameter space is rather narrow ||3q]. 

The detailed numerical analyses ]r5| , |18| , pO] , |34[|39| ] eluci- 
dating the influence of the spatiotemporal behavior of the 
concentration field on various properties of TW states, e. 
g., on the variation of flow amplitude, frequency, and 
mixing with heating rate have clearly shown that the 
success of a model description sensitively hinges upon 
the representation of the concentration field. It has to 
capture the essence of the spatiotemporal structures fol- 
lowing from the combined action of strong nonlinear ad- 
vection and weak diffusion on the one hand and the gen- 
eration of Soret induced concentration currents by tem- 
perature gradients on the other hand. A model that re- 
produces with few degrees of freedom all essentials of the 
bifurcation behavior of flow amplitude, frequency, and 
mixing is presently not available — neither in the form of 
coupled amplitude equations nor in Galerkin type form. 
The respective reasons for their deficiencies are discussed 
in the text. 

The present paper aims at filling this gap. We present 
a few-mode Galerkin model which rests upon a care- 
ful analysis 1 34 3^ of the concentration balance in liquid 
mixtures and explains among others the whole TW solu- 
tion branch from oscillatory onset up to its merging with 
the upper SOC branch and the associated changes in the 
spatiotemporal structure of the states. 



We introduce the system and formulate the theoreti- 
cal task in Sec. |l[ In Sec. fl| we construct the Galerkin 



model and give a detailed account of how the concen- 
tration field is represented. The main body of the paper 
(Sec. [V ) is dedicated to an extensive discussion of the re- 
sults. Wherever possible we provide analytic expressions 
for characteristic quantities like thresholds, bifurcation 
points, and order parameters like convective amplitude, 
frequency, heat flux, and variance of the concentration 
distribution. The SOC and TW states will be compared 
in quantitative detail with simulations. Our model also 
yields unstable SW solutions. 



II. SYSTEM 

We consider a convection cell of height d. It contains 
a binary fiuid of mean temperature T and mean concen- 
tration C of the lighter component confined between two 
perfectly heat conducting and impervious plates. This 
setup is exposed to a vertical gravitational acceleration 
g and to a vertical temperature gradient AT/d directed 
from top to bottom. The fluid has a density p which 
varies due to temperature and concentration variations 
governed by the linear thermal and solutal expansion co- 
efficients a = — and f3 — — respectively. Its 
viscosity is v, the solutal diffusivity is D, and the ther- 
mal diffusivity is k. The thermodiffusion coefficient fc^ 
quantifies the Soret coupling which describes the change 
of concentration fluctuations due to temperature gradi- 
ents. 

The vertical thermal diffusion time is used as the time 
scale d^/n oi the system and velocities are scaled by n/d. 
Temperatures are reduced by the temperature difference 
AT across the layer and concentration deviations from 



the mean concentration by ^AT. The scale for the pres- 
sure is given by 



pK 



Then, the balance equations for 



equation 

mass, momentum, heat, and concentration [p^Jl4| read 
in Oberbeck-Boussinesq approximation p0| , |l8[ 



= - 



dtu = - 



Vu 

(u- V)u- V 



P 



rf3 



dtT = 
dtC = 



R<7 (T + C) 
V • Q = - V • [uT - VT] 
V- J = -V-[uC-LV{C 



(2.1a) 



(2.1b) 
(2.1c) 
(2.1d) 



Here, the currents of heat and concentration, Q and J re- 
spectively, are introduced and T and C denote deviations 
of the temperature and concentration fields, respectively, 
from their global mean values T and C. The Dufour 
effect that provides a coupling of concentration 

gradients into the heat current Q and a change of the 
thermal diffusivity is discarded in (ETq) since it is rele- 
vant only in few binary gas mixtures ||ql and in liquids 
near the liquid-vapour critical point |2^. 
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Besides the Rayleigh number R = ^^^AT measuring 
the thermal driving of the fluid there enter three addi- 
tional numbers into the field equations (2.f a)-(2.1d): the 
Prandtl number cr = v/k, the Lewis number L — D/k, 
and the separation ratio ip — — The latter char- 

acterizes the sign and the strength of the Soret effect. 
Negative Soret coupling -0 induces concentration gradi- 
ents antiparallel to temperature gradients. In this situa- 
tion, the buoyancy induced by solutal changes in density 
is opposed to the thermal buoyancy. When the total 
bouyancy exceeds a threshold, convection sets in, typi- 
cally in the form of straight rolls for negative V'- Ignoring 
field variations along the roll axes we describe henceforth 
2D convection in an x-z plane perpendicular to the roll 
axes. 

Form and strength of convection and its influence on 
convective concentration and temperature transport are 
measured by the following order parameters: (i) The 
maximum Wmax of the vertical velocity field, (ii) The 
Nusselt number TV = (Q • e^)^; giving the lateral average 
of the vertical heat current through the system. In the 
basic state of quiescent heat conduction its value is 1 and 
larger than 1 in all convective states, (iii) The variance 



d/x,z 



(2.2) 



of the concentration field being a measure for the mixing 
in the system. The better the fluid is mixed the more 
the concentration is globally equilibrated to its mean 
value — optimally mixed, strongly convecting states 
enforce M to vanish. In the conductive reference state 
denoted by the subscript "cond" the vertical Soret in- 
d uced conce ntration gradient gives rise to a variance of 
ViC'iond)-,^ = (iv) The frequency w of a TW. 

Thus, extended TWs with a wave number k have a phase 
velocity v = uj/k. They are stationary states in a refer- 
ence frame co-moving with v relative to the laboratory 
system. 



Th e solution of the partial differential equations ( 2.1a )- 
( ^.Id ) requires boundary conditions for the flelds. We use 
realistic no slip conditions for the top and bottom plates 
at z = ±1/2, 

u(a;,z = ±l/2;t) = , 

and assume perfect heat conducting plates by 

T(a;,z = ±l/2;t) = =Fl/2 . 

Furthermore, impermeability for the concentration is 
guaranteed by 

• J = -La^ (C-V'T) (a;,z = ±l/2;t) = . (2.3) 

We should like to stress again that we restrict ourselves 
to the description of extended roll like patterns that are 
homogeneous in one lateral direction, say, y. So, we in- 
vestigate 2D states of a certain lateral periodicity length 



A = 2t: /k. In most cases we take = tt, i. e., A twice the 
thickness of the fluid layer, which is close to the critical 
wavelengths for the negative Soret couplings investigated 
here. 



III. MODE SELECTION AND GALERKIN 
MODEL 

A. Temperature and velocity fields 

The temperature field consisting of a linear conductive 
profile —2 and a convective deviation is truncated by 

r(x, z;<) = — z+To2(i)'\/2sin27rz 

+ [rii(i)e"*'=^ -l-cc] VicosTTz (3.1) 

as in the standard Lorenz model and its first exten- 
sions to convection in binary mixtures with permeable 
p6| , p6 37 1 and impermeable boundaries ]2^ , |2^ , ^ . These 
models do not provide a satisfactory representation of 
strongly nonlinear TW convection since they used a com- 
bination of concentration and temperature fields in order 
to fulfill the impermeability of the plates exactly without 
extending the temperature truncation adequately. For a 
discussion of this point sec Rcf. [Q. Here, we truncate 
the concentration itself. This approach avoids the neces- 
sity of a more complicated representation of the temper- 
ature field. 

For the velocity field we adopt an earlier successful ]3^ ] 
one-mode description 



j{x,z,t) = [wii{t)e + 



C.C. COS TTZ 



(3.2) 



Eq. ( |3.2| ) completes the Galerkin approximation of veloc- 
ity and temperature in our model. 



B. Selecting the concentration field modes 

In order to select adequate concentration modes a de- 
tailed analysis of the concentration balance and field 
structure of SOC and TW states is necessary. 



1. Lateral average of the concentration and deviation 

The first step of our analysis is to decompose the con- 
centration 

C{x, z- 1) = (C), + (C - (C),) =: Co(z; t) + c{x, z; t) 

(3.3) 

into its lateral mean Co(z;t) and the deviation c{x,z;t) 
from it. Inserting this decomposition into the balance 
equation ( p.ld| ) for the concentration and averaging it 
yields two coupled equations for Co and c 
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dtCo = -d,{wc)x + LdlC^ (3.4a) 
(St+u- V)c = a^(wc):r-w9^Co+£V2c . (3.4b) 

In both of these equations we have discarded the Soret 
couphng term Li/iV^T in the bulk of the fluid. However, 
the Soret cou phng wiU not be dropped in the bound- 
ary condition (2.3). The motivation and justification for 
this approximation is discussed in quantitative detail in 
Ref. |34[| . Here, we only mention that the basic justifi- 
cation is the smallness of the Lewis number L in liquids 
so that transport by diffusion an the Soret effect — both 
enter the balance with a weight L — are small compared 
with advection. In those regions where advection needs 
to be balanced by another transport mechanism strong 
concentration gradients are observed whereas the tem- 
perature gradient shows no such boundary layers. Thus, 
an adequate balance is assured by advection and diffu- 
sion and the additional concentration source or sink — 
the Soret effect — can be omitted in the bulk. Only in the 
impermeable boundary condition the Soret effect leads to 
a sizeable nonvanishing mean concentration gradient at 
the plates which cannot be ignored. 

In a SOC fixed point and also in a TW the lateral 
average of the concentration field is temporally constant. 
Thus, it can be calculated explicitly from (3.4a) to be 



dz' ' 



Here, the impermeable boundary condition (p.3 



a,Co(±l/2) = ?/.9,To(±I/2) 



(3.5) 



(3.6) 



relating the lateral averag es Cp and Tq has been used 
in the first integration of (3.4a) from —1/2 to z. The 
second integration is taken from to z since Co{z = 
0) = as required by the mirror glide symmetry C{x + 
X/2,z) ^ -C{x,-z) for SOC and TW states |l|,|l8|. 
When d escri bing TW as well as SOC fixed points the 
relation (3.5) can be inserted into the evolution equation 
( |3.4b| ) giving 



idt + u-V)c = LV'^c + i)Nw +{pz~Y^ 



[WC) 



(3.7) 



The solution c of this equation completely determines the 
relaxed TW and SOC concentration field whenever the 
vertical velocity field together with the Nusselt number 
is given. Just that is realized by the truncations for the 



velocity (3.2) and temperature (3.1) field in Sec. [IIA 



With c, w, and N also the TW and SOC spatial structure 
of th e lateral mean of the concentration is determined via 
(3.5). The last task is therefore to select modes for c that 



2. Symmetry decomposition and lateral mode truncation 

The mode selection of c is based among others on the 
insight gained in Ref. from a symmetry decompo- 
sition of c. The decomposition was realized p9| with 
respect to the different parities under the mirror opera- 
tions a; — > — a; and z — > — z so that the four symmetry 
classes S^^ , , S ^ , and S are obtained. The first 
(second) superscript denotes the parity under the oper- 
ation X —X (z —z). In relaxed TWs and SOCs 
the lateral coordinate x can be combined with the time 
t to X — vt so that the time derivatives in TWs can be 
replaced by —vdx and the argument t in the field c can 
be dropped. One main result of Ref. was the obser- 
vation that field components of the symmetry S are 
not needed for a quantitative description of the TW and 
SOC bifurcation topology and that the fields of the class 
are made up mainly by the zeroth lateral Fourier 
mode. Consequently, the SOC and TW concentration 
field c can be represented well by just two parts 



c{x, z) = c^^ {x, z) + c '^{x^z) 



(3.8) 



belonging to the sy mme try classes 5++ and S U sing 
the approximation ( |3.8D for c one obtains from ( |3.7D two 
equations 



V dxC + -i 
—vdxC^^ 



-w (wc++)^ = iV-i/'w + L V^c++ (3.9a) 
L 



2.- + 



(3.9b) 



for the c^'^ and c fields of fully relaxed SOCs and 
TWs. Here, the vertical velocity field w was fixe d to be- 
long to the symmetry class S^^ in the ansatz ( |3.2| ) by 
choosing the temporal phase adequately and then switch- 
ing from X to X — vt. Furthermore, we used the fact that 
the application of the advective derivatives udx and wdz 



to c 



++ 



and 



generate fields with negative vertical 



parity that do not belong to the two retained symme- 
try classes 5++ and S ^ for c. The same holds also for 



9. 



approximate the TW and SOC structure appropriately. 



The important implication of Eqs. (3.9) is that the 
lateral variation of the concentration field c{x, z) is re- 
stricted t o s in kx and cos kx if one uses the approxi- 
mations ( |3.8| ) and (3^). The reason is that w{wc^^)x 
as well as Nijjw have the lateral variation of w, i. e., 
cos kx, and no other inhomogeneity in (3.9) excites 
higher modes. We should like to stress that all these 
restrictions are based on a quantitative investigation of 
their implications. Thus, they do not endanger the suc- 
cess of our model as we will see below. 



3. Vertical variation 



The last task is to select modes for the vertical spatial 
dependence of c{x, z) in the form 
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c{x , z) = ci{z)e + c.c. 
with ci being the first lateral Fourier mode of 



(3.10) 



Nu- 



merical calculations [18 34 1 have revealed that the main 



contribution to ci(z) is made up by a part being phase 
shifted by 90° with respect to the vertical velocity field. 
With the velocity ansatz (^) and the choice that wn is 
real, this implies that Ci{z) is dominated by its imaginary 
part. This holds as long as the phase velocity v is large 
compared with the Lewis number L ~ 0(0.01). 




FIG. 1. Vertical profiles of that part of the first lateral 
Fourier mode of the concentration being in phase with the 
streamfunction (in our notation the imaginary part). Shown 
are the eigenfunctions at the Hopf bifurcation {uj — luh) and 
at the SOC-TW transition {u = 0) with arbitrary normaliza- 
tions. Dots refer to a many mode Galerkin scheme, solid lines 
to an optimal truncation by a linear combination of 1 and 
cos 27rz. Parameters are L = 0.01, a — 10, and tjj — —0.25. 

The vertical spatial dependence of Im 01(2) is investi- 
gated in Fig. [I]. Therein, we demonstrate that the "ex- 
act" solution (dotted line) obtained from a MMG scheme 
can very well be reproduced by a linear combination of 1 
and cos 2ttz, i. e., of the first two modes with the bound- 
ary condition d^c = at z = ±1/2. This holds both for 
the bifurcation out of the heat conducting state (lu = cjh) 
and for the SOC-TW transition {uj — 0). The mode am- 
plitudes of 1 and cos 27rz were optimally chosen in Fig. ^ 
in order to demonstrate the usefulness of the representa- 
tion of the first lateral Fourier mode 



Cl 



(z; t) ^ V2 [cio(t) + ci2(i) cos 27rz] 



(3.11) 



in terms of 1 and cos 27rz with two complex amplitudes 
cio(t) and Ci2(t). 



According to Eq. (3.5), the spatial dependence of the 
zeroth Fourier mode Cq{z) can be calculated by inte- 
grating the product of cos^ nz for the vertical depen- 
dence of w with 1 and cos 27rz, respectively, coming from 



c. This procedure leads to vertical modes of the form 
{2ttz + sin 27rz) and (—Anz -\- sin Airz). Their real am- 
plitudes are labelled by C02 and co4, respectively. These 
two modes C02 and co4 are constant in SOC and TW fixed 
points. To derive a model that includes temporal vari- 
ations of these modes as well, both modes are taken as 
time dependent. 

The complete Galerkin ansatz for the concentration 
field is therefore given by 



C{x, z, t) 



{l + 27r\/2[co2(t)-2co4(0]}z 
- £02(^)^2 sin 27rz + co4(t)'\/2sin47rz 



[cw{t)e~ 



ikx 



ikx 



-c.c] V2 

■c.c] \/2cos27rz . (3.12) 



In order to avoid the appearence of temperature modes in 
( |3.12D we approximate the boundary condition ( |3.6D by 
dzCo{±.\/2) — that deviates from the correct value 
by a factor equal to the Nusselt number TV = 0(1). This 
approximation can be understood as the leading term in 
an amplitude expansion of N which starts at TV = 1 . The 
exact valu e of N is of minor importance in the boundary 
condition (3.6). Only the existence of a finite slope of Cq 
at the plates is crucial. 

Similarly, a lateral variation of the vertical derivative 
of C, i. e., of c, at the plates can be seen as a higher 
order contribution that scales with the field amplitudes 
and not with 0(1). The reader can convince himself of 
the smallness of the derivatives of Imci(2) at the plates 
in Fig. g. 



C. Galerkin model 

1. Scalings 

We use r = as control parameter. Here, R'^ ~ 

I (^)^ ~ 1825.14 is the stability threshold of the qui- 
escent heat conducting state of the pure fluid with re- 
spect to disturbances of a wave number k — tt within our 
model. This is not exactly the minimum of the marginal 
curve. It is calculated as 0.9998i?° at fc = 0.98277r. But 
since we are not interested here in wave number depe- 
dencies we fix A: = tt. 

The complex amplitudes of the first lateral Fourier 
modes, / £ {wii, Tii, cio, C12}, are written in a vector 
notation 



f =(Re/ , ImfY 



(3.13) 



We scale the mode amplitudes in the following way: 

X = A^ii ' (3-14a) 
57r^ 

Y = — - — rTii , Z = — - — rTQ2 , (3.14b) 
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32V2 
Ui = — - — rcio 



25672 
Vi = — — r cq2 , V2 

lOTT 

Additionally, we introduce 



32^/2 
256\/2 



rCi2 , 



■rco4 



, 27 
cr = — cr 
14 



1 

2^ 



128 



0.6940 



(3.14c) 
(3.14d) 

(3.14e) 



2. Order parameters 

The order parameters maximal convective amplitude 
Wmax, Nusselt number A'', and mixing number M can be 
expressed by 



N 



2^^ 

/ dz(Q-e,), = 1 + — X-Y 

J-1/2 



(3.15a) 
(3.15b) 
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= 1 



75 



+ 



128 r2 

45 
64r 

1125 



|Ui| 



1 + ^ 

3 / 128 r 



1 



4 7r" 



Vo 



32768 r2 V^^ 3^^ 



75 TT^ 

16384 r2 



1^1 



(3.15c) 



Here, the Nusselt number is computed as the global spa- 
tial average of the vertical heat flux since due to the 
truncation of velocity and temperature fields in differ- 
ent bases the laterally averaged vertical heat flux which 
is conventionally used for evaluating the Nusselt number 
has a slight z-dependence. This problem occurs in all 
few mode Galerkin approximations with no slip bound- 
ary conditions, see, e. g., H,Kl[. 



3. Model 



We insert the field truncations of Se es. Ill A and 
[II B into the basic equations ( 2.11: )-( 2.1cj ) without bulk 
Sorct e ff ect and scale the mode amplitudes according to 
( 3.14a )-( 3.14d ). Then, the following model for the con- 
vection in binary fluid mixtures is obtained: 



rX = 

tZ = 
tY = 



X - Y + a0 ( Ui + 



-2{Z - X-Y) 
-Y + X{r- Z) 



(3.16a) 

(3.16b) 
(3.16c) 



tVi = 

Till = 



6L. 



-Vi 



32L, 
'15 



■1/2 + X • U 
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U2 



3 

4U2) 



-rX-lxj,-laU-lv2\X 



-rX 



5L 



10 



1, 



^2 X 



It is an extension of the standard Lorenz model 



(3.16d) 
(3.16e) 
(3.16f) 
(3.16g) 

3- The 

latter is contained in Eqs. (3.16a)-(B.16c) in a form that 
is slightly modified due to a different scaling and to the 
realistic no slip boundary conditions in our approxima- 
tion. 

This model can be looked upon as a minimal one for 
convection in binary liquid mixtures because it contains 
on the one hand the minimal description of convection in 
a pure fluid (the Lorenz model) and on the other hand 
a minimal extension for binary fluids. This extension is 
minimal since the simple extension |36 3^j2j,Q of the 
Lorenz model with only one linear mode U and one non- 
linear mode V leads to TW solutions with linear relations 
between all pairs of the three quantities (the square 
of the phase velocity), (convective intensity), and the 
Rayleigh number r. In Rcf. it has been shown that 
such pairwise linear relations result from any truncation 
of the concentration field that is limited to only one linear 
and one nonlinear mode. Clearly, these pairwise linear 
relations between v'^, w'^, and r are incompatible with the 
topology of a backwards Hopf bifurcation followed by a 
saddle node bifurcation into a branch of stable strongly 
nonlinear TWs. Thus, our incorporation of a second lin- 
ear mode U2 and its nonlinear partner mode V2 can be 
seen as a first non-trivial step in a (systematic) extension 
that goes beyond earher models p6| , ^ , p4p^ ] . 

Up to now, no few-mode model has described the bi- 
furcation topology of TWs adequately: The problem was 
not so much the backwards Hopf bifurcation but rather 
the transition to strongly nonlinear convection, the sad- 
dle node, and finally the merging of the TW solution 
branch with the upper SOC solution branch. This fail- 
ure of the earlier approximations is due to an insuffi- 
cient representation of the concentration field: It has 
not been truncated directly but rather the combination 
^ = C — TpT with the temperature field has been intro- 
duced in order to fulfill the impermeable boundary con- 
dition exactly. However, when using the combined field 

high mode representations in both, ( as well as in T 
are required as explained in [Q. By enforcing the im- 
permeability of the plates only in the lat eral a verage we 
avoid these difficulties in our truncation ( 3.12 ). 

As an aside we mention that within another minimal 
approach Knobloch and Moore |Q have deduced a model 
for free slip permeable boundary conditions. They aimed 
at a correct, analytical representation of the primary bi- 
furcation and the involved modes which is possible for 
idealized boundary conditions. However, their model 
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does not show TWs comparable with those seen in ex- 
periments and related simulations. 



IV. RESULTS 

Here, we elucidate the SOC, TW, and SW solutions of 
our model. 



A. Stationary convection 



1. Bifurcation properties of SOC states 



In the case of SOC all time derivatives in (3.16) vanish 
so that X II Y II Ui II U2 holds in the complex plane 
for the amplitudes of the laterally varying modes. With- 
out loss of generality we may therefore choose all modes 
to be real. After elimination of the temperature and con- 
centration modes the bifurcation diagram rsoc(X^) can 
be calculated as 



rsoc(X^) 



1 



1 + 



11 V 

1+X2 



1 



72 Lip 
125 X' 



^(1+X2) 



(4.1a) 



(4.1b) 



The relation (4.1) between r and X^ can be inverted, 
e. g., graphically to obtain the standard bifurcation dia- 
gram of, say, X^ vs. r. The stationary stability threshold 
fstat — ''soc(X = 0) of the quiescent heat conducting 
state follows to be 



1 



1 



^stat 



1.527| 



(4.2) 



It agrees quite well with the result Tgtat — 
1 -I- 1.538^^ of Galerkin approximations @,|l| that 
fulfill the concentration boundary condition exactly. 

As a first SOC property we can determine the type 
of the stationary bifurcation out of the quiescent heat 
conducting state. For Soret couplings smaller than 



43 a2 1 



557 r 



68 T2 
129 a 



32 r3 
215^-^ 



(4.3) 



a subcritical bifurcation is observed: dr/dX.^ < 0. The 
scaling of "0800 with L'^ agrees with earlier free slip | |2^ ] 
and no slip predictions |25| ]. For subcritical bifurcations 
the saddle node is found at 



'SOC 



. = l + —J-^^, + 0{Li;) 



1 + 1.518v/-L?/' 



(4.4a) 
(4.4b) 



being in good agreement with the numerically deter- 
mined result r|oc - 1 + 1-636V^LV? H, Eq. (4.1)]. 



Eq. (4.1b) shows that for convective amplitudes X ^ i 
the Rayleigh number corresponding to a certain ampli- 
tude square X^ deviates from that of the pure fluid 
r^=o = 1 + X^ only by terms oc Lip. This means that 
for convective amplitudes X^ ^ XItAI the bifurcation dia- 
grams of a mixture are the same as for a pure fluid. This 
equality reflects the fact that strong convective mixing in 
conjunction with diffusion equilibrates the concentration 
in the whole fluid with the exception of narrow bound- 
ary layers so that it does not influence the bifurcation 
behavior any more: The stronger the mixing the smaller 
the deviations from the pure fluid case. 



S. Fields 

In order to reconstruct the fields and with them the 
order parameters such as t he N usselt number N and the 
concentration variance M ( p.2\) we need the mode ampli- 
tudes 



Y 

Z 

V2 = 

Ui = 

Un = 



1+X2 

rX^ 
1 + X2 

1 

-^rFsoc 

-2rFsoc 
2 



■rF 



SOC 



X 

1 

X 

1 

X 

1 

X 

1 



25 

216^^ 



25 /X 

1 a — 

18 V i 



5 /X 



1 25 fx 

1 a — 

18 V L 



(4.5a) 
(4.5b) 

(4.5c) 
(4.5d) 
(4.5e) 
(4.5f) 



in the SOC fixed points. Here, we have introduced the 
quantity 

185 /X\^ 625 2 /X\*l 
1 + —a — + —a^ 



Fsoc 



185 fx 

a — 

36 \L 



648 



L 



(4.5g) 



for notational convenience. The square of M is given by 



A'L 



SOC 



1+3.905 ( — 



-f0.3040 



1+3.566 



2.224 ( 5 



0.002072 ( — 



(4.6) 



L 



0.4645 



4-1-2 



and the Nusselt number by 



SOC 



= 1 



25 X2 



18 1 



(4.7) 
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3. Comparison with numerical results 
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FIG. 2. SOC bifurcation diagrams of Nusselt number A'^ 
(a) and concentration variance M (b) vs. reduced Rayleigh 
number r. Exact (model) results are shown by dotted (solid) 
lines. 

Using these formulae we can compare the results of the 
model with " exact" ones obtained by a MMG calculation 
[Q. This is done in Fig. ||with the bifurcation diagrams 
of N and M. Since the Nusselt number [Fig. |^(a)] is 
determined by the well described temperature field the 
Nusselt number of the model deviates from the "exact" 
one maximally by 1% (at r = 1.6). The variance of the 
concentration field M in Fig. ||(b) shows that also the 
concentration field is approximated reasonably well: the 
strong mixing in stable SOCs with large velocity ampli- 
tudes (upper branch of N and lower branch of M) lead- 
ing to nearly equilibrated concentration distribution and 
nearly vanishing M is reproduced with a relative error 
of about 15%. Since with the special mode selection in 
the concentration field the model was constructed to de- 
scribe strongly nonlinear convection rather than weakly 
nonlinear states it is not surprising that M in the unsta- 
ble SOCs is reproduced only with an accuracy of about 
20% at r = 1.6. 



o 



o 




o 

II 

o 

o 



FIG. 3. Lateral (left column) and vertical (right column) 
concentration profiles of a TW with frequency lo = 2.75 (a) 
and of a SOC (b) at a Rayleigh number of r = 1.82. The 
two line types compare model (solid lines) and "exact" re- 
sults (dotted lines). Parameters are L — 0.01, a — 10, and 

= -0.25. 



Beyond these global order parameters we may also dis- 
cuss the spatial variations of the concentration field. In 
Fig. ^ we compare the SOC concentration field structure 
obatined from our model with the "exact" one from a 
MMG scheme. To that end we show vertical (right col- 
umn) and horizontal (left column) profiles. Of course, the 
model cannot describe the narrow concentration peaks in 
the lateral direction which are due to the strong bound- 
ary layer phenomena caused by the smallness of the ratio 
L/wmax = 0(0.001). Nevertheless, the model predicts 
that the concentration vanishes nearly all over the con- 
vection cell. Also in the vertical profile we see a good 
agreement when keeping in mind that the combination 
of only two modes, namely co2 and co4, can provide only 
a very rough approximation to a boundary layer. 

The quality of the approximation of the concentration 
can also be discussed by its zeroth lateral Fourier mode. 
The model predicts 



Coiz) 



t 
5 



1 y^^ ^^„^i sinn(27rz) 



-o 



L 
X 

(4.8) 



In the "exact" results Cii{z) is nearly zero in the 

bulk of the fiuid outside the boundary layers near the 
plates. It is interesting to note that an extension of the 
series in the above expression to n = cx) would yield 



1 y^^ ^y+i sin7i(27r2:) ^ 1 2ttz 



(4.9) 



for values of z e] — 0. 5, .5[ so that the linear term z in 
the square bracket of (4^) is completely cancelled. Thus, 
the exact result — Co{z) = in the limit L ^ — 
is reproduced in an optimal way, namely by giving the 
exact results for those modes that the model contains. 
A significantly improved description is possible only by 
using much more modes because the contribution from 
higher modes in (4.9) decrease only oc 1/n. 

To summarize: stationary convection in binary liquid 
mixtures is described for negative as well as for posi- 
tive ip in a semi-quantitative way by our model. Beyond 
topological details of the bifurcation diagrams even the 
peculiar spatial structures of the concentration field can 
be explained. 



B. Traveling wave convection 



1. Bifurcation and scaling properties 



For the TW fixed points of the model ( |3.16| ) with a fre- 
quency Lu we have to assume time dependences oc e*'^* of 
the complex modes X, Y, Ui, and U2 because they are 
the amplitudes of a lateral variation cx e~*'^'^. So, posi- 
tive frequencies correspond for positive wave numbers k 
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to TWs traveling to the right. The zeroth lateral modes 
Z, Vi, and V2 are time independent in TWs. We sepa- 
rate the time dependence, e*'^*, of the complex amplitude 
vectors by 



X(t),Y(t),Ui(t),U2(t) 



X,Y,Ui,U2 



(4.10) 



and use henceforth the same symbols for the time in- 
dependent prefactors. Then, by choosing the temporal 
phase, X can be taken as real while Y,Ui,U2 S C. In 
addit ion, Z, Vi , V2 € M. Inserting these solution ansatzes 
into (3.16) yields in order 0{L''^) the relations 



1 3 
— — —a— 

tp 2 a 



CT l + f72 + X2 



175 
108 ' 



1 



r!2 



1 + 



X2 



-I I 1025 „<-X\2 . 15625 „2/^X\4 
liA"'\Q.) ^ 10368" \n) 

(4.11a) 
(4.11b) 



between frequency w, amplitude X, and control parame- 
ter rxw of the TW solution. Here, 



n = 



(4.11c) 



has been introduced with t = ([3.14e ) being the in- 
trinsic time scale of the model (|3.16 ). 

Neglecting terms 0{L) in (4.11) causes the TW to 
merge at w = with the SOC solution branch of the 
pure fluid instead of with the SOC solution (4.1) of the 
mixture. Cancelling the terms of 0{L^) is allowed for all 
states with L <C or L <C X. This condition is fulfilled 
for all separation ratios away from the co-dimension two 
(CT) point where Hopf bifurcation and SOC-TW merg- 
ing fall together and none of the relations L ^ il or 
L < X can be fulfilled. The TW fixed points of our 
model can be calculated analytically without limiting to 
the order 0(L°). Since the formulae are lengthy they 
are not presented here. But they have been used for the 
calculation of the phase diagram in Fig. |^ including the 
CT point_^ 

From ( [l.llgj ) one observes first of all that TWs exist 
only for negative Soret couplings '0 < 0. This is in line 
with the absence of TWs ior ip > in numerical simu- 
lations — TWs for -0 > seen in the model o f Ref. Q 
result from low-order truncation. Relation (4.11b) al- 
lows to determine the frequency f2 of a TW with a given 
velocity amplitude X: 

rTw(X2) = i + x^ + n^{x^) 

= r^=o(X2)+02(X2) (4.12a) 
c.^(x^) ^ rTw(X2)-r^^o(X2) 
< rose - 1 ■ y ■ ) 



or 



Thus, the model predicts that the square of the frequency 
n of a TW state with a velocity amplitude X is the dis- 
tance in the control parameter r between the TW under 
consideration and the state of the pure fluid with the 
same velocity amplitude. Hence, the TW frequency is 



a direct measure of the distance of the system from the 
pure fluid, i. e., the influence of the concentration. 

Another equivalent interpretation of ( 4.11b| ) is that for 
a given fixed r the squared frequency of a TW with ve- 
locity amplitude Xtw, 



(4.13) 



is given by the difference between the pure fiuid flow in- 
tensity, X,0^Q — r—1, and the flow intensity X^^ of the 
TW in question. Thus, fl^ measures also the "vertical" 
distance in the bifurcation diagrams of X^ vs. r between 
the ■0 = pure fluid SOC solution and the TW solution 
in the mixture. 
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co(X ) / 

FIG. 4. Universal scaling relation connecting TW convec- 
tive velocity ampli tude X oc Wmax with its frequency lu. The 
model prediction (4.12b) is the identity (solid line). Numer- 
ical, finite differences (MAC) data are shown by open sym- 
bols for V e [-0.25, -0.65] and MMG data by filled symbols 
(-0 = —0.1). The symbols cover in each case the whole bifur- 
cation branch, i e., stable as well as unstable TW states. The 
inset serves as a schematic explanation of the scaling relation. 
L and a were fixed to values of 0.01 and 10, respectively. But 
the scaling relation should not depend on them as long as 
cr > 1 and L < 1. 



Eq. (4.12a) has an explicit dependence on the Soret 
coupling strength ijj since fl varies between and the 
scaled Hopf freq uency 17h. This i/j-dependence is can- 
celled by scaling ( 4.12a) wi th the Hopf frequency so that 



the left hand si de of (l4.12b ) varies for all between and 
1. Thus, Eq. (|4.12bF is a universal scaling relation for 



TW frequencies resulting from our model for small L. In 
Fig. U this prediction of the model is compared with nu- 
merical results for e [—0.65, —0.25] obtained by a finite 
difference scheme and a MMG scheme [^ . For all these 
Soret couplings the scaling relation is confirmed by the 
numerical results. Only in the case of small frequencies 
deviations are observed. They are due to the fact that 
TWs with small frequencies do not approach the SOC 
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states of the pure fluid as implied by (4.11b) but rather 
the SOC states of the mixture. These two stationary 
states differ in the Rayleigh number r by 0{L) accord- 
ing ( 4.1b| ) so that deviations of the order 0{^/L) — 0.1 
can be expected on the ordinate of Fig. ^. This devia- 
tion from the scahng relation becomes more obvious for 
weaker Soret couplings like, e. g., ip = —0.1 as shown by 
the filled triangles. These states have been computed by 
a MMG scheme . For small Soret couplings the merg- 
ing point of the SOC and TW branches is in a regime of 
small velocity amplitudes, i. e., in that part of the bifur- 
cation diagram where the differences between SOC states 
of pure and binary mixed fluids become more and more 
evident. 

It should be noticed that the scaling relation holds for 
stable as well as for unstable TWs. Furthermore, it is 
independent of the bifurcation topology, i. e., it applies 
to the form shown in the inset of Fig. ^, i. e. , for i/j > —0.4 
as well as to a topology with a bistability of slow and fast 
TWs iij < -0.4) |§. 

Bifurcation properties are best discussed by introduc- 
ing the ratio 



Wmax _ StT |X| 

V ^ 8 n 



(4.14) 



of convective and phase velocity of the TW state. Then, 
the Rayleigh number r-Tw(x^) of a TW with a velocity 
ratio X can be written in order 0{L^) as 



27 TT^ 

256 



~'~ 32 ~^ M7.X 



512' 



(4.15) 



The oscillatory stability threshold of the basic state fol- 
lows for X = as 



rTw(O) 



1 



1 



277r2 



256 1+CT 



(4.16a) 



with the numerically [g5| determined value of —3.526 • 
10-^ 

The other limit in x, namely x ^ co or w — > 0, gives 
the Rayleigh number of the SOC-TW transition: 



r* = hm rrwix ) 



1 



1+ icras^^ 

b2o 1+cr ' 



(4.18) 



However, one should keep in mind that the SOC-TW 
transition at r* with the transfer of stability from an 
SOC to a TW when reducing r is related to an instabil- 
ity of the SOC concentration boundary layer [|l9|. These 
boundary layers are caused by the smallness of the Lewis 
number L, a limit which is not systematically incorpo- 
rated in the model under consideration. The model's 
main objective is a description of strongly nonlinear TW 
convection which has for finite TW frequencies a defi- 
nite limit for small L. Thus, one should not expect a 
correct reproduction of the Lewis number dependence of 
the SOC~TW transition from a model with modes that 
do not resolve the boundary layer structure in detail. 
However, the dependence on the two other fluid param- 
eters, namely separation ratio -0 and Prandtl number a 
is given in a qualitatively correct way: strong increase of 
r* with stronger negative Soret coupling and saturation 
with i ncrea sing a as can be seen by comparing the for- 
mula ( 4.18 ) with the numerical results in |lj. Fig. 9(b) 
and 15(b)]. The SOC-TW transition points which have 
been plotted therein are also affected with a certain er- 
ror bar since the spatial resolution of the used numerical 
method was not fine enough to capture the whole bound- 
ary layer phenomena. However, their qualitative fluid 
parameter dependence has to be regarded as correct. 

The analytical form (4.15) for the TW bifurcation di- 
agram allows a simple determination of the TW saddle 
node bifurcation, namely as the minimum of rxw(x^)- K 
is given by 



' TW 



rriy (2.8687) 



0.6567 j-^V 



(4.19) 



and (4.11b) yields the Hopf frequency 

^^'''V -405.6 V 



64 



a ^ 256 'r 



27 7r2 „ 



256 



^ + 1.041 



(4.16b) 



The calculation of the CT point requires the considera- 
tion of the Lewis number dependence. If this is done one 
obtains 



CT 



1600 



459 7r2 CT 



1 



55 T 
102^ 



(4.17) 



as the separation ratio at the CT point for given Lewis- 
and Prandtl number and flxed wave number k — In 
ethanol-water mixtures with L = 0.01 and a — 10, this 
yields ilP"^ = —3.714- 10~^ being in very good agreement 



2. Fields 

The calculation of Nusselt number iV, concentration 
variance M , and concentration contrast between the two 
plates requires the computation of the temperature and 
concentration field in the TW fixed points. Their mode 
amplitudes in the TW fixed points can be expressed by 



Vi = 
Ui = 



X / L 



'i 2 / 1 2 



i] rFr 



TW 



15^ 
32' 



(4.20a) 
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Y = X(l - 
where we have introduced the quantity 



/ 1 41 2 25 4 
^--=(^+32^ +512^ 



(4.20b) 



In (4.20a) all quantities except the real parts of Ui and 
U2 are evaluated in order LP with X and taken from 



(4.11) 



3. Small amplitude expansions 

Before discussing the order parameters themselves we 
should like to show that they cannot be expanded as 
power series in the distance 



(4.21) 



from the onset of convection up to values where strongly 
nonlinear TW convection is observed. To see this, let us 



rewrite Eq. (4.21) by using (4.15) to display the relation 
between e and explicitly 



aiX^(l + Q2X^) ^ aiX^(l + a2X^) 
1 + 5ix2 + b2X^ 



1-F 



1 



(4.22) 



Here, ai^2, ^1,2, and b'l 2 are amplitude independent real 
numbers. Note in particular that the functional relation 
between control parameter e and reduced order parame- 
ter X = Wma.x/v is given by a rational function. A similar 
relation has also been found from a fit to MMG and finite 
difference numerical results 1201. The radius of conver- 



gence of a small amplitude expansion of (4.22) in powers 
of x^ is given by 



Xc 



min 1 6- 1 

1=1,2 



This quantity depends on ^ and a so that for all negative 
val ues of tp (TWs are observed only for tp < according 
to ( tl.llaD ) 



< 



328 
'25 



1481 
1681 



0.8052 



which is the absolute value of that node of -Fxw(x^) 
( 4.20b| ) with the smallest absolute value. Then, the ra- 
dius of convergence is calculated in the variable u'max/'i' 
as 



= Xc 

V J c 



0.8973 



(4.23) 



i. e., near that point in the bifurcation diagram which has 
been identified in ||20[| as the transition between weakly 
and strongly nonlinear convection. This point, namely 
Wmax — V, where areas of closed streamlines first occur, 
has also been identified as the radius of convergence for a 
small amplitude power series expansion of different order 
parameters (see Ref. pO[|). 

Thus, our model supports the notion that experimen- 
tally observed TW convection in binary liquid mixtures 
cannot be described by weakly nonlinear models as, e. g., 
complex Ginzburg-Landau amplitude equations (GL) in- 
cluding various ad hoc quintic extensions that have been 
proposed. Being used out of their validity range they can- 
not be trusted to reproduce, e. g., the relations between 
frequency to, mixing M, flow intensity w^-^^y. or Nusselt 
number N , and the thermal driving r. Typically already 
the simpler relation between w^^^ or N and r is wrong 
on the upper TW branch — not to mention the more 
sensitive relations between w, M, and r. Also results for 
LTWs based on this approach have to be questioned: 
The spatiotemporal field properties under the envelope 
being closely related to those in extended TWs are not 
captured properly. The main drawback of these GL ap- 
proaches is the insufficient representation of the role of 
the concentration field in these strongly nonlinear states. 

A first step towards a better incorporation of the con- 
centration field into the GL framework was the intro- 
duction of a long wavelength concentration mode with 
characteristic time scale cx L by Riecke The ap- 

proximation is that L is of the same order as the distance 
from onset which, however, does not apply to all exper- 
imental LTWs. Additionally, those parameters {hi and 
/i3 in Ref. [^) that could immediately lead in extended 
TW states without large scale lateral variation to a fi- 
nite mean concentration mode (possibly at the expense 
of stabilizing terms of unphysical fifth order) have been 
dropped in the LTW calculations |^^. While this ap- 
proach in its present form does not seem to generate the 
spatiotemporal field structure of TWs under the LTW 
envelope it is a promising step forwards. Incorporation 
of impermeable boundary conditions and separation of 
diffusive [L) and critical (e) time scales and thus the in- 
corporation of an additional concentration mode seems 
necessary to guarantee the aforementioned relations be- 
tween Lo, M. N. and r. 



4- Comparison with numerical results 

For a quantitative comparison with numerical MMG 
results we present in Fig. ^ bifurcation diagrams 
of the square of convective amplitude (a), the Nusselt 
number (b), the TW frequency (d), and the concentra- 
tion variance as a function of the frequency (c). They 
show that all characteristic features of the TW bifurca- 
tion scenario are captured by the model: subcritical Hopf 
bifurcation, saddle node bifurcation, stable upper branch 
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of strongly nonlinear TWs, merging of TWs and SOCs 
at w = in the strongly nonlinear regime, and drastic 
reduction of the concentration contrast with decreasing 
frequency. As an aside we mention that a model using 
the same number of modes, but numerically determined 
ones is similarly successful 39 1. This has to be contrasted 
with earlier analytical few mode approximations for TWs 
in binary mixtures which reproduced only the backwards 
Hopf bifurcation. 
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FIG. 5. TW bifurcation diagrams of the square of the con- 
vective velocity amplitude lu^ax (a), the Nusselt number A*' 
(b), and the frequency uu (d), each vs r. In (c) the reduced 
concentration variance M is plotted vs. frequency. 

The linear and weakly nonlinear bifurcation proper- 
ties, i. e., the onset of convection and the initial slope, are 
modelled by our approximation with high accuracy. The 
characteristic bump in the "exact" bifurcation diagrams 
of Fig. ||(a), (b), and (d) that occurs at the transition 
X ~ 1 from weakly to strongly nonlinear convection is not 
reproduced by our model. The absence of this fine struc- 
ture in the bifurcation diagram of our model is due to 
neglecting Fourier modes higher than the first lateral one 
in the concentration. We have explicitly checked this by 
determining the bifurcation diagrams for fields restricted 
to their zeroth and first lateral Fourier modes but with 
full vertical resolution. This is also clear from studying 
the changes in the concentration at that point: Starting 
at the onset, only harmonic lateral variation is observed 
up to a velocity ratio x — Woiuk/v ~ 1. For larger x, 
the harmonic profile gets more and more deformed by 
the occurrence of plateaus Fig. 2]. They refiect ho- 
mogenized concentration distributions in the regions of 
closed streamlines. The description of this equilibration 
requires higher lateral Fourier modes. 

The next characteristic point is the position of the TW 
saddle node bifurcation: As a consequence of the drop- 
ping of higher lateral Fourier modes and the above dis- 
cussed implications, the saddle lies at too high frequen- 
cies or too low amplitudes, i. e., in a too weakly nonlinear 
part of the bifurcation diagram. On the upper (lower) 



branch of the N vs. r {w vs. r) curve the TW frequeny 
of our model has too high values [Fig. ^(d)] and the zero 
frequency TW end point at r* lies above the " exact" one. 
This is caused by the fact that slow TWs are boundary 
layer dominated. This feature is not fully reflected in our 
model. 

Since the concentration changes significantly with fre- 
quency it is appropriate to discuss the relation between 
concentration variance M and TW frequency uj [Fig. |^(c)] 
rather than the relation between M and r thereby elimi- 
nating partly the errors in our uj vs. r curve of Fig. |^(d). 
The prediction of the model for M{iu) agrees very well 
with the "exact" curve in Fig. ^(c). This is once more a 
hint that the concentration field is globally treated in an 
adequate manner. Additionally, the relation between M 
and w is a second universal, ^/'"independent, scaling rela- 
tion when scaling the frequency with its value at the Hopf 
bifurcation. This is done in pG, Fig. 5b] with numerical 
data. 

The spatial variation of the concentration field in a 
TW is shown in Fig. |(a). To measure the quality of 
the model we compare its results with numerically ob- 
tained fields. As described above, it is convenient to se- 
lect for this procedure two TWs with the same frequency 
but different Rayleigh numbers. The "effective" value 
of the harmonic lateral profile [left part of Fig. ||(a)] in 
the model corresponds well with the " exact" plateau-like 
concentration distribution. In the vertical profile [right 
part of Fig. ^(a)] , even a slight building up of a plateau 
can be observed. Its mean height is approximated by the 
height of the lateral profile explaining the differences in 
the heights of the vertical plateaus. The strong variation 
of the model concentration along the plates is an artefact 
of only approximately fulfilling the impermeability of the 
plates. 

In the actual TW states, the concentration at the 
plates is nearly constant so that also the contrast between 
them is nearly constant. Thus, an appropriate quantity 
to compare is the laterally averaged concentration con- 
trast at the two plates. Our model predicts that 
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depends only on the velocity ratio Wmax/v- This relation 
is checked in Fig. ^ by comparing the model prediction 
with numerical finite differences and MMG results. The 
two limits, namely the basic state with Wmax and 
the stationary state v —^ (x oo) are reproduced very 
well by our approximation. The transition between them 
takes place for % e [1, 10]. 
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FIG. 6. Concentration contrast in a TW state between the 
plates as a function of the velocity ratio x ~ Wma^^/v from 
the model (solid line) and from numerical computations (fi- 
nite differences: open circles, many mode Galerkin; filled 
lozenges). Parameters are L — 0.01, a = 10, and = —0.25. 



C. Standing wave convection 

Besides SOC and TW convection there exists a third 
type of convection that bifurcates out of the ground state: 
It is standing wave convection occurring at a Hopf bifur- 
cation simultaneously with TW convection. As long as 
linear states are considered SWs are just a linear super- 
position of a right and left traveling wave. The selection 
of the nonlinear pattern leading to either a TW or SW 
is governed by a general principle |^^: stable solutions 
exist directly at a Hopf bifurcation only if both SW and 
TW bifurcate supercritically. Then, the pattern with the 
largest initial slope of the amplitude is selected. Using 
thisprinciple and the numerically obtained initial slopes 
of pSi it can be inferred that in all liquid mixtures TWs 
are stable when bifurcating supercritically, i. e., for ijj 
larger than the tricritical value V'tw °^ —L^- Thus SWs 
cannot be observed directly in experiments with liquid 
mixtures. Their investigation requires to stabilize that 
unstable fixed point and possibly destabilizing the sta- 
tionary fixed point having the same symmetries. Nev- 
ertheless, SWs represent a generic convection pattern in 
binary fluid mixtures. They occur in particular as tran- 
sients in the evolution of a strongly nonlinear TW out of 
the unstable, supercritically heated ground state. 

Up to now, only weakly nonlinear properties as the 
initial slope of SWs have been discussed for binary fluid 
mixtures with experimentally realistic, i. e., imperme- 
able boundary conditions. Since our model ( 3.16| ) has 
described both SOCs and TWs in an adequate way we 
think it is worthwhile to equally investigate the bifurca- 
tion properties and time dependence of SWs in its frame- 



work. The computation of the SW fixed points using the 
full field equations is a problem which has not been solved 
so far. 



1. Solution procedure 

In SWs the phases of the complex modes X, Y, Ui, 
and U2 are time independent so that the time derivative 
i s par allel to the mode itself, e. g., X || X. Via the model 
( |3.16D this leads to X || Y || Ui || U2, i. e., all modes 
may be chosen real without restriction of generality. 

In the linear modes X, Y, Ui, and U2 only odd mul- 
tiples of the basic frequency w of the SW occur whereas 
in the nonlinear modes Z, Vi, and V2 only even multiples 
exist. This allows to expand the time dependence of the 
mode amplitudes in the following way: 



[X(t),Y(t),Ui(i),U2(t)] 

= J2 [x("),Y("),Ui("\U2^"^ 

ri=l,3,5,... 

[Z{t),V,{t),V2{t)] 



c.c. 



c.c. (4.25a) 



(4.25b) 



These series can be inserted in the model ( 3.16| ) , the mode 
amplitude X*^^^ can be chosen real by the arbitraryness of 
a common phase in time, and the resulting nonlinear sys- 
tem of algebraic equations can be solved in the unknown 
variables {uj, X^^) , Re X^^) , ImX^^) ,...}. 



2. Bifurcation diagram and temporal behavior 

In Fig. 1^ a bifurcation diagram of the time averaged 
Nusselt number in SWs (dotted line) is given for the stan- 
dard ethanol-water fluid parameters together with the al- 
ready discussed diagrams for SOCs (solid line) and TWs 
(dashed line). The inset shows the frequency of SW and 
TW. At the oscillatory onset of convection {rose — 1.345) 
both TWs and SWs bifurcate subcritically as it has been 
predicted by |^^. With decreasing frequency a weak 
maximum in the time averaged Nusselt number of the 
SWs is observed (r ~ 1.19) before a saddle node occurs 
at r ~ 1.12. Then, it increases with further decreasing 
frequency up to a absolute maximum of about 1.1 before 
approaching the SOC branch. The SW branch has to 
be connected with the SOC branch since a SW with fre- 
quency is a SOC. The exact merging point of SWs and 
SOCs cannot be determined by means of (11251) since with 
decreasing frequency the time dependence of the modes 
gets mo re an d more anharmonic so that in the Fourier 
ansatz ( [l.25| ) more and more modes have to be included. 
The evolution of this anharmonicity with decreasing fre- 
quency is indicated in Fig. |8| where the time dependence 
of the velocity field mode X within one period of oscilla- 
tion is plotted. 
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FIG. 7. Complete bifurcation diagram of the Nusselt num- 
ber (time averaged for SWs) of all types of convection con- 
nected with the ground by a primary bifurction: stationary 
convection (SOC), traveling (TW) and standing (SW) waves. 
The inset shows the frequency bifurcation diagrams of TWs 
and SWs. 
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FIG. 8. Time dependence of the velocity field amplitude X 
over one oscillation period in SWs with four different frequen- 
cies. Parameters are L = 0.01, a — 10, and tjj — —0.25. 

For the computation of the SWs shown in Fig. |^ and 
Fig. ^ temporal Fourier modes up to lOw have been used 
and compared with states calculated with modes up to 
20(jj. Between wh and uj = 0.4 the results of both time 
resolutions agree in Nusselt and Rayleigh number better 
than 1% for fixed frequency. However, states with uj < 
0.4 shown in Fig. M and Fig. may be erroneous. 



D. Phase diagram 



-0.01 -0.001 




FIG. 9. Phase diagram in mixtures with L = 0.01 and 
a = 10 in a double logarithmical plot r — 1 vs. tp. SOC 
properties are shown by dotted lines (saddle node t^qq: open 
triangles up, bifurcation rgtat: open triangles down). Solid 
lines correspond to TWs (Hopf bifurcation rose : filled squares, 
SOC-TW-transition r*: open lozenges, saddle node r^^: 
filled circles). Stable TWs exist within the shaded area. 

The dependence of the bifurcation topology of SOCs 
and TWs on the separation ratio is summarized in the 
phase diagram of Fig. |^ which can be compared directly 
with the "exact" results described in Q Fig. 7]. The 
characteristic points of the SOC bifurcation diagrams, 
namely the stationary onset rgtat and the saddle node 
^soc agree very well with the "exact" results. The same 
holds without restriction for the oscillatory onset of con- 
vection rose and the TW saddle node bifurcation r^yj- 
The '(/'"dependence of the SOC-TW transition point r*, 
i. e., the strong increase with increasing the negative cou- 
pling strength, is only reproduced qualitatively: r*{'ip) 
runs too flat in the interval [—0.3, —0.01] so that the point 
where r* = rose is given as V' — —0.014 instead of —0.14. 
This is due to the neglection of the influence of boundary 
layers in the concentration field. The same is true for the 
merging of r* with the TW saddle node, i. e., that separa- 
tion ratio beyond which no stable upper TW branch can 
be observed: ip ~ —0.001 instead of —0.008. Neverthe- 
less, the model predicts a -^-interval where the TW sad- 
dle node can be seen below the SOC saddle node, as it has 
been discussed in the framework of the "exact" results 
in ||, Fig. 5d]. Furthermore, the tricritical SOC-TW 
transition Fig. 5e] at which the TW branch merges 
vertically with the SOC branch, is predicted to take place 
on the unstable SOC branch, too. For more negative ■0 a 
TW saddle appears. Then, stabilization of TWs at that 
saddle and subsequent destabilization, probably towards 
modulated TWs, occurs before the now unstable TWs 
end on the unstable lower SOC branch jsj. Fig. 6]. 
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V. CONCLUSION 

We have investigated roll like 2D convection in bi- 
nary liquid mixtures with negative Soret coupling. Then, 
three types of extended convective states occur which are 
connected via a primary bifurcation to the quiescent heat 
conducting state: stationary convection (SOC), traveling 
waves (TW), and standing waves (SW). One objective 
of our paper was to derive a model describing the com- 
bined SOC-TW bifurcation topology and the character- 
istic spatiotemporal behavior of the concentration seen 
in numerical simulations and experiments. The most im- 
portant TW bifurcation features are (i) a backwards Hopf 
bifurcation, (ii) a saddle node giving rise to stable TWs, 
and (iii) the merging of the TW solution branch with 
the SOC branch. Along the TW branch, (iv) both the 
phase velocity and the variance of the concentration field 
decrease monotically in the same way. 

Model — To derive such a model we started with an 
approximation for velocity and temperature fields simi- 
lar to that in the standard Lorenz model |^^. We used 
however a more realistic improved version with a no slip 
velocity field. To select and motivate an appropriate 
representation of the concentration field we relied on a 
systematic analysis of the concentration balance equa- 
tions: the structure of the field components occurring in 
a symmetry decomposition was investigated , a sepa- 
ration into lateral mean fields and deviations thereof was 
used, and the effect of the Soret coupling in the bulk 
of the fluid layer and at the plates were quantitatively 
assessed for liquid mixture parameters jsj] . The concen- 
tration field truncation derived from these investigations 
consists of two linear and two nonlinear modes. The re- 
sulting Galerkin model can be looked upon as minimal 
for the description of convection in binary liquid mix- 
tures since it contains the Lorenz model, i. e., the sim- 
plest truncation for the pure fluid, and a minimal exten- 
sion for mixtures. Previous extensions [ p6p7|j2^ , ^ of 
the Lorenz model were too simple and therefore failed 
in reproducing strongly nonlinear properties like, e. g., 
the TW saddle node. The present model is the first an- 
alytically manageable approximation showing the above 
stated four characteristic properties (i) - (iv). 

Stationary convection — Good agreement in the prop- 
erties of the SOC branch was found: stability threshold, 
fluid parameters at the tricritical bifurcation, and posi- 
tion of the saddle node bifurcation. Furthermore, the 
approach of the SOC branch of a binary mixture to that 
of a pure fluid could be discussed in the limit of large 
convective amplitudes. Even the spatial variations of the 
concentration field which is boundary layer dominated in 
SOCs are reproduced in a way allowing good quantita- 
tive agreement in the concentration variance. This holds 
for both negative and positive ip. 

Traveling waves — The main results have to be seen 
in the description of TWs. Here, the model shows that 
TWs occur only for negative Soret couplings in agree- 



ment with all numerical simulations and experiments. It 
predicts that the TW frequency is a direct measure for 
the "distance" of the system from the pure fluid, i. e., 
for the influence of the concentration field. The "dis- 
tance" can directly be read off the bifurcation diagrams 
of flow intensity versus Rayleigh number in two equiv- 
alent geometric ways. This insight yielded a universal 
scaling relation between phase velocity, convective veloc- 
ity, and degree of mixing of a TW which was confirmed 
in an impressive manner by different numerical methods 
analysing the full field equations. The derived scaling 
relation holds for all TW states, stable or unstable inde- 
pendent of the bifurcation topology |Q . 

Linear convective properties are reproduced by the 
model with high accuracy: oscillatory stability thresh- 
old, Hopf frequency, and CT point. The same holds for 
the Rayleigh number at the saddle node bifurcation. Our 
model shows also a SOC-TW transition to SOCs at the 
upper stationary stable branch and its dependence on 
the separation ratio and Prandtl number is qualitatively 
correct. However, its Lewis number dependence is un- 
physical since the concentration boundary layers which 
are responsible for the SOC-TW transition ||l^ are rep- 
resented in the model only in an incomplete way. Oth- 
erwise, the spatial structure of the concentration field in 
TWs — frequency dependence of the concentration con- 
trast between the two plates and building up of plateaus 
(in the vertical direction) — is modelled in a quantita- 
tively correct way by our truncation. 

The model allows to pinpoint the breakdown of an ex- 
pansion of the TW solution as a power series in the dis- 
tance from the onset of convection up to values where 
stable, strongly nonlinear TWs are observed in experi- 
ments. The related radius of convergence of the model is 
close to the numerically determined one pO[ | marking the 
transition from weakly to strongly nonlinear states [ pO[ . 
Hence, complex Ginzburg-Landau equations should not 
be expected to yield reliable quantitative results for lo- 
calized and extended TW states. 

The results of the stationary and traveling states are 
brought together in a phase diagram whose good agree- 
ment with numerical simulations can directly be in- 
ferred. Only in the SOC~TW transition remarkable de- 
viations are observed. 

Standing waves — Our model gives insights into non- 
linear SW solutions in binary mixtures. It confirms ear- 
lier weakly nonlinear results |^,^ like the initail slope. 
In addition, it becomes possible to follow the SW branch, 
which is everywhere unstable, up to regions with strongly 
nonlinear oscillating amplitudes. A numerical determina- 
tion of the unstable SW solution of the full field equations 
is still lacking and an observation of these states, say, by 
a control process is an experimental challenge. 
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